How continuous quantum measurements in finite dimension are actually discrete 



O 

o 



> : 

00 . 

\o ■ 
o ■ 

(N ■ 
O 

o : 

Oh: 

a ■ 



X 



Giulio ChiribellEEl and Giacomo Mauro D'Ariancjll 
QUIT Group, Dipartimento di Fisica "A. Volta" and INFM, via Bassi 6, 27100 Pavia, Ital^ 

Dirk Schlingemanrll 

Institut fiir Mathematische Physik, Technische Universitdt Braunschweig, 
Mendelssohnstrafie 3, 38106 Braunschweig, Germany 
(Dated: February 1, 2008) 

We show that in finite dimension a quantum measurement with continuous set of outcomes is 
always equivalent to a continuous random choice of measurements with only finite outcomes. 

PACS numbers: 03.65.Ta, 03.67.-a 



When we measure the spin component along a mag- 
netic field with a Stern-Gerlach apparatus, for spin 1/2 
particles we have only two possible outcomes: spin up 
and spin down. This measurement is perfectly repeat- 
able, and can perfectly discriminate between the two 
orthogonal states |t) and |J,). It is possible, however, 
to design an experiment with more than two outcomes, 
which discriminates optimally — though not perfectly — 
among three or more non orthogonal states. Indeed, 
a four-outcome measurement on a two-level system is 
needed in the eavesdropping of a BB84 cryptographic 
communication [l[, or in the lab to perform an informa- 
tionally complete measurement^, which determines the 
quantum state from the measurement statistics. 

What about performing a measurement with a con- 
tinuous set of outcomes? This is the case of a mea- 
surement designed to optimally determine the "di- 
rection" of a spinj^, similarly to what we do in 
classical mechanics. Such a measurement produces 
a probability p{n)dn of the spin direction falling 
within the solid angle d n around the direction n = 
(sin 9 cos (f), sin 9 sin cos 9) . Indeed the measurement 
of direction must be feasible^ — though in-principle 
inaccurate — otherwise Quantum Mechanics would fail in 
describing what we normally observe in the macroscopic 
world. Actually, this is not the only interesting ex- 
ample of continuous-outcome measurement on a finite- 
level system: in fact, measurements of this kind have 
an endless number of applications, e.g. optimal state 
estimation^, optimal alignment of directions^ and ref- 
erence frames optimal phase estimation Q, and opti- 
mal design of atomic clocks Q. 

In this Letter we establish a fundamental property of 
quantum measurements with continuous set of outcomes, 
namely that for finite level systems any such measure- 
ment is equivalent to a continuous random choice of mea- 
surements with only finite outcomes. This means that 
any physical quantity measured on a finite dimensional 
system is intrinsically discrete, while the continuum is 
pure classical randomness. For a spin 1/2 particle, this 
fact is well illustrated by the simple observation that the 
optimal measurement of direction can be equivalently re- 



alized by a customary Stern-Gerlach experiment where 
the magnetic field is randomly oriented. We emphasize 
that, in general, the discretization of physical quanti- 
ties does not involve just von Neumann observables, but, 
more generally, finite measurements with a number of 
outcomes larger than the Hilbert space dimension. Mor- 
ever, using the main result we show that any continu- 
ous measurement that optimizes some convex figure of 
merit (e.g. maximizing the mutual information or the 
Fisher information, or, alternatively, minimizing a Bayes 
cost[l, 01) can be always replaced by a single measure- 
ment with finite outcomes, without affecting optimality. 

Let us start by briefly reviewing the general theoreti- 
cal description of measurements in Quantum Mechanics. 
Consider a quantum system (with Hilbert space TL of di- 
mension dim(7Y) = d < oo), which undergoes a measure- 
ment with random outcome u, distributed in the outcome 
space fi. The probability distribution of the outcomes for 
any possible state p of the system depends on the spe- 
cific measuring apparatus used, and is given by a positive 
operator-valued measure (acronym POVM), namely the 
probability that the outcome falls in the subset B C is 
given by the Born rule p{B) = Tr[pP(i3)], where P{B) 
is a nonnegative operator with normalization condition 
P{fl) = I Uj. In the special case in which the measure- 
ment is finite, a random result i from a set of possible 
outcomes {i — 1,2,..., N} is returned with probability 
Pi — Tr[pPi], Pi ^ being nonnegative operators with 
normalization condition X^iLi Pi — ^■ 

Before presenting the main result, in order to help in- 
tuition, we briefly analyse two simple prototypes of con- 
tinuous measurement: the optimal measurement of the 
"spin direction" for a spin 1/2 particle, and the optimal 
measurement of a phase shift. 

The measurement of direction for spin 1/2 particles is 
given by the POVM Q 



PiB) = / ^\n){n\ 
Jb 27r 



(1) 



where \n) is the eigenvector of n- J with eigenvalue -t-1/2, 
J being the spin operator. It is simple to see that this 
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measurement is equivalent to the randomization 

P{B)^ [ ^E^-Hb) , (2) 
JS2 47r 

where dn/(47r) is the uniform probabihty distribution 
over the unit sphere and iJ^") is the POVM 

E'^^\B)^XB{n)\n){n\+XB{-n)\-n){-n\ . (3) 

[xsin) is the characteristic function of the set B: 
XB{n) = 1 for n G B, xs(") = otherwise.] The POVM 
represents a measurement of direction based on a 
Stern-Gerlach setup with magnetic field oriented along n: 
if the apparatus outputs "up" , one assigns to the spin the 
direction n, if 'down" , one assigns —n. With this data- 
processing, the probability of observing the spin within 
the region B is nonzero only if B contains at least one 
of the directions ±n. The continuous POVM ([T]) is then 
interpreted as a Stern-Gerlach measurement performed 
with a direction n of the magnetic field randomly chosen 
in the unit sphere. 

Another example of continuous measurement is that of 
phase estimation, where one wants to measure the phase 
shift (f> e [0, 27r) experienced by a quantum state under 
the action of the unitary evolution — exp(iAf0), with 
N ~ J2t=o {!")} orthonormal basis for H. The 

optimal POVM is given by [i| 

and is equivalent to the randomization: P{B) = 
^l"' ^E^'t'\B), where E^'l''^ is the POVM E'^'t'^B) = 

2 EtJo Xb(0„ + 0) |0„ + +^\,K = ^. 

We will now show that all continuous measurements in 
finite dimension can be always interpreted in an analo- 
gous way, namely as a continuous random choice of mea- 
surements with finite number of outcomes. More pre- 
cisely, we will prove the following 

Theorem 1 For any POVM P{B) the following decom- 
position holds 

P{B) = / dxp{x) E^^'^iB), (5) 
Jx 

where x £ X is a suitable random variable, p{x) a prob- 
ability density, and, for every value of x, E^^^ denotes a 
POVM with finite support, i.e. of the form 

E{B) ^Y.^b{lo,)P, (6) 

i=l 

{uji S VL] being a set of points, and {Pi} being a finite 
POVM with at most d^ outcomes fT^ I . 



A POVM E of the form of Eq. ^ is in turn nothing 
but the continuous data-processing of the finite POVM 
Pi, with function of the outcomes /(i) = uJi'. if the ap- 
paratus outputs i, then one assigns to the measurement 
the outcome uii. The decomposition ^ shows that the 
continuous POVM P is achieved by randomly choosing 
a classical parameter x G X and then performing the 
finitely-supported POVM E^^\ depending on x through 
the finite POVM {P/"'''} and through the points {uj'l'^^}. 
Operationally, this corresponds to the following recipe: 
i) randomly draw a value of x according to p{x); ii) 
depending on x, measure the finite POVM p(^), thus 
getting the outcome i\ Hi) for outcome i, assign to the 
continuous measurement the outcome As a first 

consequence, this simple recipe shows that, contrarily to 
a rather common belief (see, e.g. jl3|). continuous quan- 
tum measurements in finite dimension are as feasible as 
the discrete ones. 

The decomposition of the measurement of the " spin di- 
rection" given by Eq. ([2]) provides a concrete example of 
decomposition ([5]). In particular, the finitely-supported 
POVM in Eq. © is illustrated in Fig. |T]for n = k. 
Notice that, in general, there may be different random- 
ization schemes yielding the same continuous POVM: as 
an example. Fig. [1] illustrates another finitely-supported 
POVM that allows one to reproduce the measurement of 
direction by simply randomizing the orientation of the 
Cartesian axes. 




FIG. 1: Left: Illustration of the POVM S*"' in Eq. © as 
an example of finitely-supported POVM S'^' in Eq. ((6]). In 
this specific example the outcome space Q, is the unit sphere 
n = S^, the dimension of the Hilbert space is d = 2, and 
only two out of the four terms in Eq. ((6]) are nonvanishing, 
namely Pi = |fe)(fe|, P2 = \ - k){~k\, P3 ^ P4 = 0. The 
probability of finding the spin direction in a region i? C 
is zero for R missing the two poles, as B in the figure, and 
is possibly nonzero for R as A. Right: Another example 
of finitely-supported POVM for d = 2, corresponding to a 
SIC (symmetric informationally complete) POVM 14]. The 
POVM is made of four elements Pi corresponding to Ui at 
the vertexes of a tetrahedron. The probability can be nonva- 
nishing only if the region R contains at least one of these four 
points, such as in ^, whereas it is always zero in situations as 
in B. 

We now derive the main result. We fix both the quan- 
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turn system and the outcome space ft, and consider the 
set T of all possible POVMs for these. This is a convex 
set, since given any two POVMs P' and P", their con- 
vex combination P(^) = XP' + (1 - A)P" for A G [0, 1] 
is still a POVM, namely the whole segment joining P' 
and P" is contained in T. The extremal points of the 
convex set 7 are those POVMs that cannot be written 
as convex combination of two different POVMs. Stated 
differently, a POVM P £ CP is not extremal if and only 
if it is the midpoint of a segment completely contained 
in y, i.e. if and only if there exist two distinct points 
P\P" ey, P' P" such that P = \{P' + P"). This 
is equivalent to the existence of a direction Q ^ Q and 
a positive number e > such that P + tQ S ? for any 
t G [— e, e]. The standard name for the direction Q in 
convex analysis is perturbation. Here the perturbation Q 
is a function that associates to any subset P C $1 an op- 
erator Q{B), fulfilling the three requirements: i) Q{B) 
is Hermitian for any subset B C fl; ii) Q{fl) = 0; Hi) 
P{B) + tQ{B) > for any P G and for any t G [-e, e]. 

If there exists a nonzero perturbation Q for P, then 
P is non extremal: using this criterion, we now establish 
that the extremal POVMs must necessarily have finite 
support, namely they must be of the form of Eq. ([6]). 

The proof takes advantage of the following: 

Lemma 1 Every POVM P G T admits a density with 
unit trace, namely for any POVM P there exists a finite 
measure jiiduj) over fl such that 

P{B) = f M(d^) M(a;) , (7) 

with M{ll}) ^ and Tt[M{uj)] — 1 fj,-almost everywhere. 

Proof. Consider the finite measure /i(dtj) defined by 
^(P) = Tr[P(P)], yB C n . Since P(P) > 0, one has 
P(P) ^ Tr[P(P)] / = n{B) I, namely P(P) is domi- 
nated by the measure fi{B). This implies that P admits 
a density M{uj) with respect to fi{du!). Clearly, the den- 
sity M{uj) has to be nonnegative /x— almost everywhere. 
Moreover, for any B C fl one has Jgii{duj) = /x(P) = 
Tr[P(P)] = JB^^{dLu) Tt[M{uj)], whence Tr[M(w)] = 1 
/i-a.e. ■ 

Thanks to this Lemma, we can represent any POVM 
P G IP using its density M{u!) as in Eq. ([7]). To prove 
that an extremal POVM must be of the form ([6]) it is 
enough to show that for extremal POVMs the measure 
/i(dw) must be concentrated on a finite set of outcomes 
{uJi, . . . jUJd^}, i.e. fJ.{B) = for any set P C not con- 
taining anyone of the points Ui. We recall the definition 
of the support of a measure ^.{Auj) as the set of all points 
w G 57 such that /^(P) > for any open set P containing 

LJ. 

Lemma 2 Let P G 7 be a POVM and fi{diu) the mea- 
sure defined by fJ.{B) = Tr[P(P)]. If P is extremal, then 



the support of fi^duj) is finite and contains no more than 
d^ points. 

Proof. Suppose that the support contains more than 
points. In this case, one can take d^ -{- 1 points cj^ G 51 
in the support and d^ -\- 1 disjoint open sets UiC ft, 
i — 1, . . . , + 1, such that tOi G Ui for any i[15[. As 
a consequence, the space L°°{fl,fi) of integrable func- 
tions f{Lij) that are bounded /i— almost everywhere has 
dimension at least d^ + 1 (indeed, the characteristic func- 
tions xUi i^) are a set of -I- 1 bounded and linearly 
independent functions). Then, consider the matrix ele- 
ments fijiij) = {i\M{uj)\j}, where M{uj) is the POVM 
density of Eq. d?]), and \i), \ j) are elements of an or- 
thonormal basis for 7i. Since the operators M{oj) are 
nonnegative with unit trace a.e., the functions fij{i-o) 
are bounded a.e., namely fij G L°°{ft,ii) yi,j. More- 
over, since the space L°°{rt, /i) has dimension larger than 
d^, it must contain at least one function gito) that 
is linearly independent from the set {fij}. Using the 
Gram-Schmidt orthogonalization procedure, such a func- 
tion g can be always chosen to be orthogonal to all fij, 
namely J^^{duj) g* {to) fij (uj) = Vi,j. Finally, since 
f*j{uj) — fji{(^) yi,j, such a g can be also chosen to be 
real. Now, we claim that the Hermitian operators Q{B) 
defined by 

Q(P) = / fi{duj) g{uj) M{uj) (8) 
Jb 

provide a perturbation for the POVM P. Indeed, we have 
Q{il) = as the matrix elements {i\Q{fl)\j) are zero for 
any i,j: 

{i\Qmj} = / f,{dcj) gicj)miu;)\j) (9) 

= f fi{dLo)g{Lo)f,jiLo) = . (10) 

Moreover, since g G L°°{il,fi), there exists a positive 
number c such that |3(w)| ^ c < oo a.e., thus implying 
that the operators M{lu) [1 + tg{Lu)] are a.e. nonnega- 
tive for any t G [— e, e], with e — l/(2c). Hence, inte- 
grating over any subset P, we obtain that the operators 
P(P) -\-tQ{B) are nonnegative, namely Q is a perturba- 
tion. Finally, Q is nonzero, otherwise taking the trace 
of Eq. ([8]), and using that Tr[M((jj)] = 1 a.e. we would 
get = Tr[g(P)] = /sM(dw) g{uj) VP, thus implying 
5 = 0, which is not possible by definition of g. In con- 
clusion, if the support of /i(dw) contains more than d^ 
points, then the POVM P has a nonzero perturbation, 
whence it is not extremal. I 

Lemma [2] establishes that an extremal POVM has nec- 
essarily the form of Eq. ^ , namely it can be realized by 
measuring a finite POVM Pi and declaring measurement 
outcome uJi. Using this fact, we readily obtain the proof 
of the main Theorem: 
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Proof of Theorem [TJ Due to the standard Krein- 
Milman theorem of convex analysis, any point of a com- 
pact convex set is a continuous convex combination of 
points that are either extremal or limit of extremals. On 
the other hand, it is simple to prove that the set J" of all 
POVMs is compact [l^, and that any limit of extremal 
POVMs is still a POVM of the form For a de- 

tailed mathematical proof see Ref. • ■ 

We now want to explore some consequences of decom- 
position ^ for optimization of POVM's and for Quan- 
tum Tomography. Optimizing a quantum measurement 
consists in finding the POVM P that maximizes the value 
of a figure of merit T[P] — e. g. mutual or Fisher in- 
formation, average fidelity, or any Bayes gain. In all 
these cases J^[P] is convex, i. e. J^[AP' + (1 - X)P"] < 
XJ'iP'] + (1 - \)T[P"] for any A G [0, 1]. Suppose now 
that a continuous POVM P is optimal for !F. Combining 
convexity of T with Eq. ([5]) one has 



Keyl, L. Cattaneo, and S. Albeverio for many useful dis- 
cussions in the preliminary stage of this work. 



^max = ^[P] < / dxp{x) T[E<^-^] < 
J X 



(11) 



which implies T\Ef^^'^\ — ^max for any x except at 
most a set of zero measure. This means that all the 
finite POVMs S(^) are equally optimal: in particular, 
/or any optimal continuous measurement there is always 
an optimal measurement with finite (no more than cfi) 
number of outcomes. In special situations some ex- 
plicit algorithms to find optimal finite measurements are 
known [lil [Til. [2o|. In particular, Ref. [3] shows that in 
many cases the minimal number of outcomes is larger 
than d = dim(7i). Combined with the above result, 
this fact definitely proves that the quantum discretiza- 
tion cannot rely solely on von Neumann measurements. 

Regarding Quantum Tomography, using the present 
analysis we can make mathematically precise the com- 
mon intuition that an informationally-complete contin- 
uous measurement is equivalent to a Tomography scan 
made of a random choice of observables — or more gen- 
erally POVMs. Indeed one can estimate the ensemble 
average of any operator A by using the two data pro- 
cessing /yi(w) and fA{i^i^^) for continuous POVM and 
tomography, respectively, as follows 



A 



X 



Axp{x)y jA{u:r')P, 



{x)-.-p(x) 



In conlusion, in this Letter we showed that continu- 
ous quantum measurements can be always realized by 
performing finite measurements depending on a random 
classical parameter. Physical properties, such as spatial 
orientation and time, are then intrinsically discrete when 
measured on finite level quantum systems. 
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